Quantum Dynamics of Lorentzian Spacetime Foam 
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ABSTRACT 

A simple spacetime wormhole, which evolves classically from zero throat ra- 



cr 

S-h ! dius to a maximum value and recontracts, can be regarded as one possible mode of 



fluctuation in the microscopic "spacetime foam" first suggested by Wheeler. The 
dynamics of a particularly simple version of such a wormhole can be reduced to that 
of a single quantity, its throat radius; this wormhole thus provides a "minisuper- 
space model" for a structure in Lorentzian-signature foam. The classical equation 
of motion for the wormhole throat is obtained from the Einstein field equations and 
a suitable equation of state for the matter at the throat. Analysis of the quantum 
behavior of the hole then proceeds from an action corresponding to that equation of 
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motion. The action obtained simply by calculating the scalar curvature of the hole 
spacetime yields a model with features like those of the relativistic free particle. In 
particular the Hamiltonian is nonlocal, and for the wormhole cannot even be given 
as a differential operator in closed form. Nonetheless the general solution of the 
Schrodinger equation for wormhole wave functions, i.e., the wave- function propa- 
gator, can be expressed as a path integral. Too complicated to perform exactly, 
this can yet be evaluated via a WKB approximation. The result indicates that 
the wormhole, classically stable, is quantum-mechanically unstable: A Feynman- 
Kac decomposition of the WKB propagator yields no spectrum of bound states. 
Though an initially localized wormhole wave function may oscillate for many classi- 
cal expansion/recontraction periods, it must eventually leak to large radius values. 
The possibility of such a mode unstable against growth, combined with the observed 
absence of macroscopic wormholes, suggests that stability considerations may place 
constraints on the nature or even the existence of Planck-scale foamlike structure. 
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I. INTRODUCTION 



A quantum description of gravitation is one of the most eagerly sought goals 
of present-day physics. Approaches to such a description may be loosely classified 
as "grand schemes" and "small schemes." The former are attempts at a compre- 
hensive quantum theory of gravitation, such as supergravity or superstring theory, 
canonical quantization of general relativity with or without new variables, et cetera. 
The latter are attempts to analyze the quantum-gravitational physics of particular, 
simple systems, in hopes of understanding features expected to emerge from a com- 
plete, general theory. These small schemes include the extensive "minisuperspace" 
quantum-cosmology program [1] , the numerous calculations of the effects of Planck- 
scale spacetime wormholes and "baby universes" on the constants of macroscopic 
physics [2], and the present work, in which we seek to probe the dynamics of such 
a Planck-scale "foam" structure itself. An outline of this investigation has been 
published previously [3]; here we present the calculations and results in detail. 

The concept of "spacetime foam," first suggested by Wheeler [4] some 35 years 
ago, has become standard lore in our quantum picture of gravitation: On scales 
of the Planck length, quantum fluctuations of the spacetime geometry become so 
dominant that spacetime takes on all manner of nontrivial topological structure, 
such as wormholes and handles; smooth, simply connected spacetime only emerges 
as a classical limit on larger scales. Yet in the intervening decades our understanding 
of quantum gravity still has not advanced sufficiently to prove or disprove this 
famous conjecture. 
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The difficulties inherent in such a proof, or in any detailed treatment of space- 
time foam, are well known. Topological change in a Lorentzian manifold necessarily 
entails the complications of singularities or degeneracy of the geometry, or of closed 
timelike curves [5] . To circumvent these problems one may consider Euclidean man- 
ifolds instead [1] . Euclidean quantum gravity is physically distinct from Lorentzian, 
unlike ordinary quantum field theory in which the two formulations are equivalent 
via analytic continuation. Indeed on a quantum level it might be that spacetime is 
fundamentally Euclidean, with familiar Lorentzian spacetime only a classical limit 
attained far from the Planck regime. But Euclidean quantum gravity has its own 
formidable problems [6], including the failure of the Euclidean action to be posi- 
tive definite or bounded below, the problem of interpretation, and the matter of 
recovering Lorentzian spacetime. 

Both versions of spacetime foam continue to be subjects of great interest. 
Much recent work on Euclidean foam has explored its possible role in determining 
the fundamental constants of nature [2]. The suggestion that Lorentzian wormholes, 
extracted from microscopic foam and suitably enlarged, might conceivably be used 
for space or time travel [7] has inspired a great deal of recent and ongoing work 
on issues of causality and consistency in spacetime physics. Both sets of ideas 
have profound significance for fundamental physics, but both regard the spacetime 
foam itself as a given, neither probing the questions of its existence or structure 
substantially. 
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This work addresses the dynamics of Lorentzian spacetime foam, and the 
implications of that dynamics for its existence and nature. We begin by envisioning 
Lorentzian spacetime, filled on Planck scales with all manner of wormholes and other 
structures — modes of topological fluctuation — continually winking into existence, 
persisting for microscopic time periods, and pinching off. The moments of creation 
and disappearance of these structures, i.e., the actual points of topological change, 
may or may not require a Euclideanized treatment; our analysis does not attempt 
to deal with these points. 

Classical wormhole geometries can be used to model such a picture. In general 
such wormholes, persisting for arbitrary lengths of time (in particular, sufficient to 
allow passage of matter or energy), must involve stress-energy distributions which 
violate the weak energy condition — negative energy density in some reference frames 
must appear somewhere in the wormholes' throats. Whether this suffices to rule 
out wormholes on macroscopic scales is unknown at present, but it is not expected 
to be a problem for microscopic structure, where quantum effects could readily give 
rise to such stress-energy [7,8]. Moreover it does not appear to present a barrier 
to the growth of wormholes from Planck scales to large sizes, as the present work 
indicates. 

To make the analysis tractable we employ a brutally simplified wormhole 
model to represent one mode of spacetime-foam fluctuation. A class of worm- 
holes can be constructed by excising the world tube of some surface from a 3 + 1- 
dimensional spacetime region and joining this to another such region, with like ex- 
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cision, at that surface [9] . The join, the throat of the resulting wormhole, contains 
a surface layer of stress-energy specified by the Einstein field equations expressed as 
junction conditions [10]. Here we take the external regions of such a wormhole to 
be flat, empty Minkowski space, the throat a sphere of time- varying radius [9]. Our 
model is thus an extreme version of a wormhole in which the spacetime curvature 
in the throat is much greater than that in the regions surrounding the mouths. Our 
simplifications reduce the geometric degrees of freedom of the spacetime to just 
one, the throat radius. Thus our analysis is a minisuperspace model of a wormhole, 
analogous to the minisuperspace models used in quantum cosmology [1]. In fact 
our model is simplified even further than those: We treat the matter providing 
the stress-energy on the wormhole throat via an equation of state, rather than as 
a separate dynamical field. The quantum/gravitational dynamics of the complete 
wormhole system becomes then the quantum mechanics of one variable, the throat 
radius. The result we seek is the evolution, in time as defined in the wormhole's 
flat exterior regions, of a wave function for that radius. 

In the absence of a general theory of quantum gravity, there is no prescription 
from first principles for quantizing a system like this wormhole. The existence 
of an external time, in terms of which the evolution of the system is naturally 
to be described, makes unsuitable the usual time-independent Wheeler-de Witt 
equation [11]. Instead we choose to "solve the constraint" explicitly, i.e., to impose 
the Hamiltonian constraint classically, to reduce the phase space of the system. 
Combining the classical junction conditions at the wormhole throat — consisting of 
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the Hamiltonian constraint plus a dynamical equation — with an equation of state for 
the matter on the throat yields an equation of motion for the throat radius alone. We 
then construct an action for the dynamics in the reduced phase space, i.e., depending 
only on the radius and its time derivatives, for which that equation of motion is 
the Euler-Lagrange equation. This can be done in a variety of ways, amenable 
to various forms of quantization. One such action, associated with gravitational 
dynamics in a straightforward way, is obtained from the scalar curvature of the 
wormhole spacetime. The wormhole quantum mechanics described by this action 
is treated in detail in this paper. Alternative actions and their implications for 
wormhole dynamics are examined in a subsequent paper [12]. 

The wormhole quantum mechanics which follows from the "curvature action" 
is in some respects similar to the Newton- Wigner first-quantized description of a 
relativistic free particle [13,14]. For both the Hamiltonian is non-polynomial in the 
momentum, i.e., is a nonlocal operator in the coordinate representation. Indeed, for 
the wormhole the Hamiltonian is a rational function of "velocity," which is related to 
the canonical momentum via a transcendental equation which cannot be inverted 
in closed form. This unconvential form of the Hamiltonian makes it impossible 
to deduce the behavior of the system merely by examining the "potential" term; 
explicit solution for the evolution of the wave function is needed. Moreover the 
corresponding Schrodinger equation for wormhole wave functions cannot be written 
explicitly as a differential equation. It can, however, be solved, at least formally: 
The wave-function propagator is obtained by using the action in a Feynman path 
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integral. Exact evaluation of this integral is problematic; the action is complicated, 
and questions of the appropriate measure and skeletonization (factor-ordering prob- 
lems) , as well as the issue of the correct class of paths over which to integrate, arise. 
But the integral can be evaluated in a WKB approximation, as a sum of contri- 
butions from classical paths and small fluctuations about them. The measure and 
factor-ordering problems are thus avoided. It is crucial, though, to include all the 
appropriate paths. The analogy to the relativistic particle suggests that spacelike 
as well as timelike and lightlike paths must be included in the integral [14,15]. Also, 
since negative throat radii are undefined, a suitable boundary condition must be 
imposed on the wave functions, i.e., on the propagator, at zero radius. Zero radius 
acts as a reflecting barrier, giving rise to classical paths — extrema of the action— 
which "bounce" there one or more times. (Such a boundary condition incorporates 
our neglect of any topological change at zero radius.) The WKB version of the 
propagator, then, consists of a sum of several terms, contributions of paths with 
different numbers of bounces; cancellation of terms effects the boundary condition 
at zero, just as for a particle in a half space. Explicit expressions for all the terms 
are obtained here. With these the evolution of a wave function is calculated via a 
convolution integral, evaluated numerically. 

The results suggest that even with a matter equation of state chosen so that 
the classical evolution of the throat radius is bounded, the wormhole is quantum- 
mechanically unstable. This is not unprecedented — classically stable black holes, 
for example, are unstable to Hawking evaporation. In this case a Feynman-Kac 
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decomposition of the wormhole propagator yields no spectrum of bound states. 
The calculated evolution of an initially localized wave function follows a classical 
trajectory for many expansion/recollapse/bounce cycles, but the wave function must 
eventually leak to arbitrarily large throat-radius values. Its behavior thus resembles 
that for a particle initially confined to a well with finite, though perhaps high and 
wide, walls, or an a particle in a heavy nucleus. This simple wormhole geometry 
appears to represent a mode of spacetime-foam structure unstable against growth 
to macroscopic size. 

Though subject to the limitations of our model, viz., the simplifications and 
approximations we have applied, and the uncertainties associated with our quan- 
tization scheme, such a conclusion is potentially very significant. It indicates that 
Planck-scale gravitational physics might be constrained by observations on much 
larger scales: The apparent absence of any topological structure to spacetime on 
scales accessible to observations, ranging from particle physics to astronomy, could 
imply either that not all modes of topological fluctuation are possible or that mi- 
croscopic (Lorentzian) spacetime foam does not exist altogether. 

The wormhole model and its associated classical and quantum dynamics are 
presented in Sec. II below. Evaluation of the wormhole propagator in the WKB 
approximation is shown in detail in Sec. III. The implications of the result for 
wormhole wave-function evolution are discussed in Sec. IV. Conclusions and caveats 
follow in Sec. V. Units with G = c = K = 1 are used throughout. Sign conventions 
and general notation follow those of Misner, Thorne, and Wheeler [16]. 

9 



II. MECHANICS OF THE WORMHOLE MODEL 
A. Classical dynamics 

A spherically symmetric "Minkowski wormhole" [9] provides a very simple 
model of a mode of topological fluctuation in Lorentzian-signature spacetime foam. 
The classical geometry of such a wormhole is constructed by removing a ball of 
time-dependent radius r = R(t) from two Minkowski-space regions and identifying 
the two boundary world-tubes, making them the throat of the wormhole. The 
Einstein field equations are satisfied trivially in the flat, empty exterior regions. At 
the throat those equations take the form of the junction conditions [10] 

S) = ±-{K)-8)K™\, (2.1) 

where Sj is the stress-energy tensor of a surface layer on the throat, and the right- 
hand side is the jump discontinuity in the throat extrinsic curvature Kj, minus its 
trace K™. Here the nontrivial components of these equations are 

1 1 

Stt = ~2^R (l- £2)1/2 (2 ' 2a) 



and 



i / R R 2 R \ 



See = to U-#)V2 + (i - R 2 ?' 2 ) ' (2 ' 2b) 

overdots denoting derivatives with respect to Minkowski-coordinate time t (in a 

frame "centered on the hole," i.e., in which spherical symmetry is maintained). The 

throat coordinates here are proper time r and polar angle 9; proper and coordinate 

time are related via dr = (1 — R 2 ) 1 ^ 2 ^. Once a surface-layer equation of state, 

relating surface density a = S TT and pressure p = Sgg/R 2 , is specified, Eqs. (2.2) 
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give the classical equation of motion for the wormhole. As is to be expected [7,8], 
Eq. (2.2a) indicates that negative energy densities appear. 

For simplicity we specify the equation of state as a property of the model, 
rather than derive it from some more fundamental description of the surface layer. 
It could be chosen to make the equation of motion extremely simple. The choice 
p = —a/2, for example, would imply R = 0. The dynamics of such a model is trivial. 
Classically the wormhole throat either sits at fixed radius or expands or contracts 
linearly. Its quantum behavior is the same as that of a free particle; a wave function 
initially concentrated about some R value will disperse to infinity. But this is not a 
model suitable for describing structures in spacetime foam, fluctuating into and out 
of existence. Rather we seek an equation of state such that the equation of motion 
describes expansion from zero radius to some maximum value and recollapse. The 
simple "dust" equation of state, p = 0, implies the equation of motion RR—R 2 + 1 = 
0; the solutions of this are sine functions, with the desired behavior. This equation, 
however, is not suited to the action construction we use below. (It can be treated 
by other methods [12].) Instead we choose 



p = -a/4 , 



(2.3) 



which yields the equation of motion 



2RR - R 2 + 1 = . 



(2.4) 



The corresponding classical trajectories are parabolic: 




(2.5) 
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where a and to are constants. These solutions likewise behave as desired. 

Thus by restricting consideration to Minkowski- wormhole geometries, and 
combining the classical initial-value (Hamiltonian) constraint (2.2a) with the dy- 
namical equation (2.2b) and the equation of state (2.3), we reduce the infinity of 
gravitational and matter degrees of freedom of the wormhole spacetime to one, the 
throat radius R. This is thus an extreme form of minisuperspace model. 

The equation of motion (2.4) can be obtained as the Euler-Lagrange equation 
for extrema of the action 

1 + R 



S = J \RRhi 



2R) dt . 



(2.6) 



1-R 

This choice of reduced-phase-space action is motivated by the fact that it coincides 
with the Einstein-Hilbert action for the wormhole spacetime. The wormhole's scalar 
curvature TZ is nonzero only on the throat: 



n = 2{K^\5{p) 
4 



R 2 {1 - R 2 ) 1 / 2 _ 





1 + R 


RR In 




1-R 



2 dt \ 



1 + R 


)] 


1-R 





S(P) , 



(2.7) 



where p is r—R(t) in one exterior region and R(t) — r in the other. The corresponding 
gravitational action — the integral of 1Z plus a surface term [17] eliminating the total- 
derivative term — is just Eq. (2.6). 

The Hamiltonian dynamics of the wormhole, on its reduced phase space, fol- 
lows immediately from the action. With Lagrangian L given by the integrand in 
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Eq. (2.7), the canonical momentum is 




1 + R 2RR 
1-R 1-R 2 



) 



(2.8) 



The Hamiltonian is 



H = PR — L = 



2R 



(2.9) 



1-R 2 



The Hamiltonian cannot be expressed in terms of R and P in closed form, since the 



It should be noted that the "energy" corresponding to this Hamiltonian is not 
to be identified with the ADM mass of the wormhole, which is zero by construction. 
To thus fix the "energy" would reduce the equation of motion to first order, and the 
classical initial-parameter space to one dimension, making canonical quantization 
impossible. In fact the Hamiltonian (2.9) corresponds to the classically conserved 
"mass" 2(2nR?/ 2 a) 2 appropriate to the equation of state (2.3). [The Einstein field 
equations imply the conservation law d(AnR 2 a)/dt + pd(4irR 2 ) /dt = 0; Eq. (2.3) 
then implies that _R 3 / 2 cr, hence H, are classically conserved.] This can take on 
various values for different wormhole states. In this respect our treatment contrasts 
with the analyses of Hajfcek et. al. [18] of the quantum mechanics of a collapsing 
dust shell. Those authors take the Hamiltonian to correspond to the ADM mass of 
the system, which may take on different values in different states. They treat the 
classically conserved proper mass of the shell as a fixed parameter. 

Of course the wormhole action is not uniquely defined by the equation of 
motion specifiying its extrema. A variety of actions, of other than Einstein-Hilbert 



transcendental relation (2.9) between R and P/R cannot be inverted explicitly. 
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form, can be constructed corresponding to Eq. (2.4) — or to equations obtained 
from other equations of state. Actions can be obtained which avoid some of the 
complexities encountered with the geometric action (2.6), e.g., the complicated 
kinetic term, and the transcendental relation (2.8). The construction of such actions 
and their implications for wormhole quantum mechanics are examined in another 
paper [12]. The rest of this paper treats the quantum dynamics obtained from the 
geometric action. 



B. Quantum dynamics 

In this minisuperspace model, then, the quantum wormhole is to be described 
by a wave function ip(R,t). With the Hamiltonian constraint (2.2a) imposed clas- 
sically, i.e., incorporated into the reduction of the wormhole's phase space, no 
Wheeler-de Witt equation imposes "time independence" of ip [11]; instead its evo- 
lution in wormhole-exterior time t is generated by the Hamiltonian (2.9). 

There is no hope of treating the corresponding Schrodinger equation directly: 
The transcendental nature of relation (2.8) and the form of H imply that the Hamil- 
tonian contains all powers of the momentum P, or as an operator on ip, of d/dR. 
It is therefore a nonlocal operator, similar to that for a relativistic particle, but 
rather more unwieldy [14]. Operator ordering in H is problematic as well. But the 
solution of the Schrodinger equation can be treated. The general solution, i.e., the 
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time evolution of any wave function, can be given in terms of a propagator: 



with C denoting the class of paths over which the integral is defined. 

The appropriate class of paths is an important issue. The example of the rel- 
ativistic particle [14], and more general quantum-gravitational considerations [15], 
suggest that all paths linking the initial and final points should be included— 



to "negative energy" contributions, in terms of the Hamiltonian (2.9).] Including 
spacelike paths implies that G need not vanish outside the "light cone" \R — Rq\ — t, 
i.e., that it may admit "acausal propagation." 

Moreover, the wormhole geometry is only defined for nonnegative R, hence 
only paths with R(t) > should be included. This restriction can be implemented 
as for a particle confined to a half space, i.e., as if there were an infinite poten- 
tial wall at R = 0. This implies the boundary condition i/j(0,t) = 0, hence that 
G[R, t; Ro 7 0] should vanish for R or Rq zero. (More precisely, ip must entail no cur- 
rent in the —R direction at R = 0, so boundary conditions more general than this 
are possible. For simplicity we do not treat such here.) By imposing this condition 
we exclude consideration of topology-changing processes, i.e., wormhole creation or 




(2.10) 



The propagator G is given formally by a Feynman path integral 




(2.11) 



spacelike as well as timelike and null. [Spacelike paths, with \R\ > 1, correspond 
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disappearance, at R = 0. The treatment of such processes might require a "second 
quantized" framework (actually "third quantized," outside the restrictions of a min- 
isuperspace model), rather than the "first quantized" (actually "second quantized") 
description used here. Their bearing on our results is discussed further in Sec. V. 
The condition on ip and G gives R = here the role of a reflecting boundary for 
wormhole wave functions. 

III. EVALUATION OF THE PROPAGATOR: WKB APPROXIMATION 

Exact evaluation of the path integral (2.11) is out of reach: The integral, with 
action (2.6), is non-Gaussian; the operator-ordering problems involved in canoni- 
cal quantization with the Hamiltonian (2.9) would reappear as ambiguities in the 
skeletonization of the path integral; the appropriate measure on the space of paths 
is unknown. But over a suitable range of time intervals and throat-radius values it 
should be possible to calculate the propagator via a WKB approximation. Then the 
path integral is taken to be dominated by the contributions of classical trajectories — 
local extrema of the action — and nearby paths. The approximate propagator takes 
the form 



The exponential factors arise from the classical paths R c \ between the initial and 
final values. The prefactors come from Gaussian integrals over small deviations 



dent of the choice of path-integral measure. Specifying the regime in which this 




(3.1) 



Paths 



from these, giving semiclassical corrections [19]. To that order G?( WKB ) is indepen- 
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approximation is accurate is problematic, since we have neither an exact solution 
for comparison nor even an explicit form of the Schrodinger equation. This question 
is considered further in Sec. V below. 



The first step in evaluating G^ WKB \ then, is finding the appropriate classi- 
cal paths. More than one type of path always contributes. There is exactly one 
parabolic trajectory of form (2.5), and taking only nonnegative radius values, be- 
tween any initial and final points (i?o, 0) and (R, t). Specifically, this trajectory has 
parameters 

V\ 1 ' 2 Ro + R 



a 



t 2 



( t 2 \' 



(3.2a) 



and 



to 



( t 2 \ 1/2 
t [R R+-\ -R 



RqR + 



t 2 \ 1/2 R + R 



(3.2b) 



and action 



S (Q) = _^L 
2a 



2 

l-^(t 2 -3tt + 3^) 



R\ 
+ T ln 



2-a(t- to) 



"T ln 



2 + at 



2 — ato 



(3.3) 



2 + a(t-t ) 

Because paths are restricted to the half-space of positive radius values, there are 
also classical trajectories which are piecewise of form (2.5), but which "bounce" one 
or more times at zero radius; extremizing the action determines the bounce times. 
Such a path which bounces just once, at time t\, has action 



= - j[ti(ti + 2Ro) + (t- h)(t -h + 2R)} 



+ \ [r 2 In 



R 
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R 



-R 2 In 



h-Ro 



Rq 



(3.4) 



For the bounced path to contribute to the WKB propagator (3.1), this should be 
extremal with respect to variation of t\. That condition is 



(t - ti) 2 (ti - Rq) - t\{t - h - R) = . 



(3.5) 



Each real solution of this equation with t\ G [0,£] corresponds to a contributing 
trajectory. Since the left-hand side of Eq. (3.5) changes sign between t\ = and t\ = 
t there is always one such solution, and may be three. A classical path bouncing 
more than once, at times t\, . . . , t n , has action 



S (n) = - k[h(ti + 2R ) + J2(t 3 - tj-i) 2 + (t- t n )(t -t n + 2R)\ 

1=2 



+ \ [ R 2 In 



R 



t — tn — R 



-Rl In 



ti-Ro 



Ro 



(3.6) 



That this be extremal with respect to all bounce times is equivalent to the conditions 



(n - l)tj - (h - R )(t n - h) = , 



j — 1 , 

tj =t\-\ -{t n - t\ t 

n - 1 



2 < j < n - 1 



and 



(n -l)(t- t n ) 2 - (t n - h)(t -t n -R) = 



(3.7a) 



(3.7b) 



(3.7c) 



Any real solution to the coupled quadratic equations (3.7a) and (3.7c) with < 
ti < t n < t yields an n-bounce path contributing to c( WKB ). For a given n there 
may be zero, two, or four such. Conditions (3.5) and (3.7a-c) are equivalent to 
requiring that the parameters a for all the parabolic segments of a bouncing path 
be equal, which is to say that the wormhole "bounces elastically." 
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For any set of initial and final coordinates (Rq,0) and (R,t) there is a maxi- 
mum number of bounces n max , beyond which there are no appropriate solutions to 
Eqs. (3.7a) and (3.7c) and no contributions to g( wkb ). For example, in the case 
Ro = R, the four solutions of Eqs. (3.7a) and (3.7c) are 

t[ a) = (t- \jt 2 - ARnt] /(2n) 

1 /(2n) , 



M - 



(2n - l)t - \A 2 - 4Rnt 



(3.8a) 



tf } = (t+ \jt 2 - ARnt^j /(2n) 
= [(2n - l)t + ^t 2 - 4i?nt 



/(2n) , 



(3.8b) 



and 



(c) 



1 

M 



= (t + 2R + ^/t 2 -ARnt + AR 2 ^j /[2(n + 1)] 



(2n + l)t - 2R - \jt 2 - ARnt + AR 2 /[2(n+l 
(t + 2R- \A 2 - ARnt + AR 2 ^j /[2(n + 1)] 



/[2(n+l)] 



(3.8c) 



(3.8d) 



(2n + 1)* -2R+ \jt 2 -ARnt + AR 2 
For n between 2 and t/(AR), if any, all four solutions are real and in the desired 
interval; four n-bounce paths contribute to G^ WKB \ For any n between t/(AR) 
and (t 2 + AR 2 )/(ARt) the (a) and (b) solutions are complex, and there are no cor- 
responding classical paths, but the (c) and (d) solutions characterize two n-bounce 
paths contributing to the propagator. For n greater than (t 2 + AR 2 )/(ARt) there 
are no n-bounce classical paths. (The values Ro = R = 0, for any £, are a limiting 
case, in which four n-bounce paths contribute for all n greater than 1.) It must 
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be emphasized that these particular bounds for n, and the classification of the so- 
lutions (a), (6), (c), and (d), cannot be applied for arbitrary unequal Rq and R, 
but there are always some bounds on the n values of contributing paths. This is 
because between bounces the paths are of form (2.5), for which there is a fixed ratio 
between maximum radius and the interval between zero crossings. For any t value, 
then, there is time for only so many bounces high enough to reach a given Rq or 
R value. 

The evaluation of G^ WKB ^ proceeds from these results. Expression (3.1) even- 
tually takes the form 

q(wkb) [Rt t . ^ Q]= £ £ p (n) ( £LJM exp , (3 . 9) 

with n the number of bounces, as above, and k labelling the n-bounce paths, 
the range of k depending of course on n. The classical actions Sj^ are given 

(n) 

by Eqs. (3.3), (3.4), and (3.6). The functional determinants are: 

F<°> = - 32 2 /( " t} , , (3.10a) 
4 + a 2 t (t-t ) v ' 



with a and to from Eqs. (3.2); 



F (i) _ hJih - R o) 2 

t 2 + 2R t - 6tti + Qt\ + 2(R - R )t! 



(3.10b) 



with t\ from Eq. (3.5); 



F {n) _ tjjtn -h)/(ti - Rq) 



t[{4n 2 -2n- l)h - (2n - l)t n - 2(n - 1)R ] + R[t n - (2n - l)*i] ' 
+2n{t\ + t 2 ) - \n 2 t x t n + R [(2n - l)t n - t{\ 

(3.10c) 
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for n > 2, with t\ and t n from Eqs. (3.7a) and (3.7c). The phase factors are 
pj^ = +1 if has the same sign as and pj^ = — 1 if it has the opposite 
sign. This choice of phases ensures the boundary conditions G( WKB )[0, t; R , 0] = 
and G( WKB )[.R, t; 0, 0] = 0, hence the condition -0(0, t) = for any wormhole wave 
function, via the cancellation of contributions from paths differing by one bounce. 
For example, as R approaches zero, the contribution of the direct (zero-bounce) 
classical path becomes equal in magnitude and opposite in sign to that of a one- 
bounce path 1 with bounce time approaching t. This phase prescription accords 
with that used, for example, for the simple problem of a particle confined to a half 
space [20]. 

Form (3.9) and the accompanying expressions show that even the approximate 

propagator for this simplified system is a complicated quantity. Some of its features 

are illustrated in Fig. 1. The propagator has singularities on the outgoing and 

ingoing or "bounced" light cones R = Ro+t, R = Ro—t, and R = t—Ro : respectively. 

(It has other singular points as well, described below.) As expected, given the 

inclusion of spacelike trajectories in the path integral as discussed in Sec. II, the 

propagator is nonvanishing outside the light cone. Thus it admits acausal, i.e., 
1 The absolute value of the argument of the logarithm is taken in the action (2.6) 

and subsequent expressions to ensure that this cancellation still occurs when the 

direct path is just timelike and the one-bounce path just spacelike, i.e., for Rq t 

and R approaching zero. 
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superluminal, propagation. Overlying these gross features the propagator exhibits 
a great deal of high-frequency structure. 

IV. QUANTUM EVOLUTION OF THE WORMHOLE 

In principle the propagator contains a complete description of the quantum 
dynamics of the wormhole. Analysis of the approximate form C?( WKB ), carried out 
in several stages, indicates that the wormhole is quantum-mechanically unstable 
against growth to large size. 

Straightforward numerical evaluation of the propagation integral (2.10), using 
propagator G?( WKB ) a s given by Eq. (3.9) and a simple initial wave function, reveals 
behavior of the wormhole over short and intermediate time scales. Results of such 
calculations are displayed in Fig. 2. For this example the initial wave function is 
a localized "wave packet," ipo = (2/-7T) 1 / 4 exp[— (Rq — 10) 2 ], all quantities in Planck 
units. For the purpose of numerical integration this is taken to be zero outside the 
interval 6 < Rq < 14. Although all terms in Eq. (3.9) have been given analyti- 
cally, the numerical integration is nontrivial due to the existence of singular points 
in G( WKB ). The poles in the propagator on the outgoing, incoming, and bounced 
light cones are treated numerically by replacing the divergent term in g( wkb ) ? in 
a small interval (denoted A) in Rq about the light cone, by its average over that 
interval. This is determined from an expansion in Rq of the appropriate term about 
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the light cone. In particular we obtain 

R-t+A / F (0)\V2 
R-t-A Urf) -pK } )^ = 

f * 



4i 



' t(2i?-t) \ 1/2 _ 
J t(2R-t) 



Am' 



x exp< 



#A 



3i 



- ^t(2i2-t)| +0(A 2 lnA) 

(4.1a) 



on the outgoing light cone, 



R+t-A \ Am 



x exp - 



t(2^ + t)lnf t(2 f A +t) U(i? + t) 2 ln / R 



RA 



R + t 



3z 



t(2i? + t)| + 0(A 2 lnA) 

(4.1b) 



on the ingoing light cone, and 
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' (t - R)A ' 
R 2 + (t-R) 2 / 



3z 



[# 2 + (t-#) 2 ]j+0(A 2 lnA) 

(4.1c) 



on the bounced light cone (where only one single-bounce, and no multiple-bounce, 
classical path contributes). In the indicated intervals, then, these singular terms 
in the propagator are replaced by 1/ (2 A) times the right-hand sides of these ex- 
pressions. We have checked that although the integrals (4.1) depend on A, the 
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final integrals (2.10) are independent of A over a range of values from 10~ 4 to 10 _1 
Planck lengths. The WKB propagator has additional singularities, caustics, where 
new bouncing paths appear, i.e., where Eq. (3.5) or Eqs. (3.7a and c) have coincident 
roots. These are avoided by evaluating Eq. (2.10) at t values which are transcen- 
dental to machine accuracy; rational values of R and Ro are taken so the t values 
at which caustics would occur are algebraic or rational. Figure 2, then, shows the 
resulting wave function tjj(R, t). Its most prominent feature is that the wave packet 
follows an essentially classical bouncing trajectory. Some rapid spreading of the 
packet appears early on, but this behavior does not seem to recur at later times. 
Numerical difficulties hinder the reliable evaluation of ip at t values larger than 
those shown: Accurate numerical integration becomes problematic as the propaga- 
tor develops ever higher- frequency behavior, and the problem of accurately solving 
Eqs. (3.5) and (3.7a and c) with large coefficients makes even calculating the prop- 
agator difficult. To determine the long-term behavior and stability of the wormhole 
a different sort of analysis is required. 

The late-time behavior of the propagator can be examined analytically in a 
special, but very useful, case. For arbitrary Rq, R, and t the explicit form of (7( WKB ) 
is intractable — the general expressions for the bounce times t\ and t n appearing in 
the and are too unwieldy even for computer manipulation! But for Rq = 
R the solutions (3.8) yield manageable results, and nearby Ro and R values can 
be handled perturbatively. The terms in the sum (3.9) are then labeled a, 6, c, 
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and d, corresponding to those solutions. In the limit t 3> -R, -Ro, with R = Ro + e 
and lei <C -Ro, the functional determinants and classical actions in those terms are: 



C(n) _ _1_ 
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for n less than (t 2 + ABjq) / {ARtf) . The n = term is a c term, given by Eqs. (4. 2d) 
and (4.2e), and the n = 1 terms are a, 6, and c terms given by Eqs. (4.2a-e). 
The order-e terms in the actions are obtained from the relation dS c \/dR = P(t) 
and Eq. (2.8) for the momentum P, evaluated on the appropriate classical path. 
Terms of order e and higher in the functional determinants give rise to subdominant 
contributions and are not needed here. 

The above expansions in inverse powers of t do not converge for n = t/(ARo). 
Indeed, the a, 6, and c terms in the propagator are singular for these parameter 
values, at which the roots (3.8a-c) coincide. But the contribution of these singular 
terms to the evolution of a wave function can still be estimated. At fixed Rq, 
such a singularity or "caustic point" occurs at intervals in t of ARq, the period 
of a classical bouncing trajectory with maximum radius value Rq. The singular 
contribution is thus associated with the classically bouncing peak seen in Fig. 2. 
Unlike the similar singularity in a harmonic-oscillator propagator [21], however, this 
singularity in (^( WKB ) does not correspond to a 5 function. Rather, in the immediate 
vicinity of the caustic, i.e., for e small compared to t~ 2 (in Planck units), the c term 
has functional determinant 

'c<">~^ + - (4.3a) 

and action 

S (n)^_ J R t_^! e 2/8 + ... _ (43b) 

The contribution to a wave function from this vicinity is therefore of order t~ 3 / 2 , 

not unity, as a 5 function would yield. However, the form of the action (4.3b) 
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suggests that a significant contribution from this term may arise from a region of 
width t -1 / 2 in e. Such a contribution might be of order unity, as Fig. 2 suggests, 
but the expansions used in Eqs. (4.3) do not suffice to calculate it precisely. The 
a and b terms do not in fact exist for nonzero e about the caustic point; they give 
no contribution here. Another set of caustic points occurs where the c and d roots 
coincide. Here again the a and b terms do not contribute, those roots being complex, 
while the c and d terms are similar in form to those given by Eqs. (4.3). 

Hence it is easiest to examine the behavior of g( wkb ) for parameter values 
such that no caustics occur at integer n values. Then the expansions (4.2) give 
approximate values for all terms. The actions vary slowest with e, i.e., with i?o, 
for the largest-n terms, so it is these which determine the width of the region in 
integral (2.10) which gives a significant contribution to a wave function. That width 
corresponds roughly to one oscillation of the slowest-varying terms. For example, 
the integral of G( WKB ) [R, t;R-e,0], with terms given by Eqs. (4.2), over the interval 
e G [—tt/R, +tt/R], is shown in Fig. 3 for fixed R and various t values. These results 
suggest that a wave function at fixed radius value continues to oscillate even for 
quite large t values, as might be expected from a sum of terms with rapidly varying 
phases. Thus they still give no clear picture of the very-long-term behavior of the 
wormhole. 

But telling features of that behavior can be extracted from the propagator as 
approximated by Eqs. (4.2). Did the Hamiltonian (2.9) describe a bound system, its 
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ground-state energy (or, for a free system, the bottom of its continuous spectrum) 
would be given by the Feynman-Kac [22] limit 

E Q = - lim - \nG[R,-ir;R,0] (4.4) 

of the propagator in imaginary time. This can be evaluated via Eqs. (4.2). The 
propagator is still a sum of terms with rapidly varying phases with respect to 
increasing r, owing to the t 2 terms in the actions. Consequently it oscillates forever; 
the limit on the right side of Eq. (4.4) does not exist. This is illustrated in Fig. 4. 
This means that the quantum wormhole possesses no spectrum of bound states. 
The expectation value of R cannot be confined in the late-time limit. As for a 
system with an inverted potential diverging to negative infinity, or a metastable 
system such as a particle confined by finite walls, the wormhole wave function must 
eventually run or "leak" to arbitrarily large radius values — the wormhole is unstable 
against eventual growth to large size. The previous results suggest that it is the slow 
leaking behavior which characterizes these wormholes. They appear reminiscent of 
a particles in a nucleus, oscillating perhaps millions of times before escaping to 
infinity. 

V. IMPLICATIONS AND LIMITATIONS 

Spherically symmetric Minkowksi wormholes [9] provide a simple model of a 

mode of topological fluctuation in Lorentzian spacetime foam, a mode apparently 

unstable against growth to macroscopic size. The quantum-gravitational dynamics 

of these wormholes is reduced to the quantum mechanics of a single variable, the 
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throat radius, by describing the matter at the wormhole throat with a suitable 
equation of state and imposing the Hamiltonian constraint at the classical level to 
reduce the phase space of the system. The corresponding reduced action is used in a 
Feynman path integral to obtain the propagator for wormhole wave functions; this 
is evaluated in the WKB approximation. The result indicates that these wormholes 
have no bound quantum states. Though their throat radii are classically bounded, 
i.e., they are classically stable, they will nonetheless grow to large size by quantum 
"diffusion." 

Many systems exhibit similar behavior. For a particle with the familiar 
quadratic kinetic term in the action, the form of the potential determines whether 
such diffusion or spreading occurs: A potential well with walls or barriers which fall 
off at large distances will allow a classically bound particle to leak out via quan- 
tum tunneling (as in the case of a decay), while one which increases monotonically 
with distance will not. For these wormholes, with more complicated action (2.6), 
so simple an analysis is not possible. The more involved examination of the worm- 
holes' quantum dynamics described here is needed to see that spreading of the wave 
function to large radii will take place. 

Our result suggests that dynamics, and stability considerations in particular, 
may be of great importance in understanding the quantum nature of gravitation. 
A definitive demonstration of the existence of an unstable mode of fluctuation in 
spacetime foam would have profound implications. Since a macroscopic structure of 
wormholes is not observed, i.e., spacetime appears to be smooth and topologically 
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trivial on all scales accessible to laboratory physics, it would imply the existence 
of a mechanism for suppressing such a mode, or even the absence of (Lorentzian- 
signature) spacetime foam altogether. Of course the present work is far from such a 
definitve demonstration; it serves to point up lines along which these matters should 
be studied further. 

The most fundamental limitation of our calculation is the restriction of the 
gravitational degrees of freedom to those of the spherically symmetric Minkowski 
wormhole, i.e., the use of a "minisuperspace model" for topological structure. In 
fact our model is even more restricted than the usual minisuperspace models [1], 
since the matter in the hole is treated not as a dynamical field but by the use of 
an equation of state. Moreover we use the particular equation of state (2.3), to 
simplify the calculations; other possible choices are considered in Ref. 12. 

We analyze our constrained model by quantization in the reduced phase space, 
as described above. In the absence of a general framework for quantum-gravity 
calculations, this method seems best suited to the problem. It does differ markedly, 
though, from the Wheeler-de Witt approach [11]. 

Furthermore, we use the particular reduced action (2.6). Other forms corre- 
sponding to the classical equation (2.4) are possible; the effect of this choice on the 
results will be examined elsewhere [12]. 

Moreover our choice of action implies a choice of Hamiltonian fundamentally 
different from that used in similar calculations [18], corresponding here to a classi- 
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cally conserved "proper mass" rather than the ADM mass of the wormhole. This 
choice raises some interesting questions. For example, since the Hamiltonian (2.9) 
is time independent, its expectation value is conserved. Its form appears to suggest, 
then, that the expectation value of the throat radius R should remain bounded. It 
might be expected, however, that any initially localized wave function must contain 
both positive- and negative- "energy" (\R\ > 1) components; this is a well-known 
feature of, e.g., the quantum mechanics of a relativistic particle [13,14]. The inter- 
play of these components could account for the eventual spreading to large radius 
values. More detailed calculations of the wormholes' quantum dynamics should 
clarify this. 

We employ the WKB approximation to evaluate the wormhole propagator. 
WKB calculations of quantum instabilities in classically stable systems, such as 
tunnelling, are well known. Here it is more difficult to be precise about the accuracy 
of the approximation. It is certainly to be expected to be valid in the late-time 
limit t R, Rq, in which the instability is manifest. But the accuracy of the 
numerical evolution of a wave function, as shown in Fig. 2, is harder to establish. 
Lacking any exact solution for comparison, we have tested the accuracy of the 
calculation via the composition relation 

G ( WKB )[#,t ; # ,0] = J G ( - WK ^[R,t;R h t 1 ]G ( - WK ^[R 1 ,t l ;R ,0]dR 1 , (5.1) 

with < t\ < t. It can be shown analytically that this should hold if the WKB 

approximation is strictly valid at the intermediate time t\. We found, however, that 

for some t\ values relation (5.1) is not well satisfied. This may indicate inaccuracy 
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of the WKB approximation, or may be due to numerical difficulties associated with 
the integration, given the rapidly varying phase of the propagator. 

A final limitation, of fundamental significance, is our implementation of the 
restriction that throat radii are nonnegative. Here we do this as for a particle in a 
half space, leading to the boundary condition ip(0,t) = 0. Other implementations 
might be used, the most general condition being only that the wave function ip 
entail no current in the —R direction at R = 0. Our condition eliminates from 
consideration any processes such as wormhole creation or disconnection at R = 
0. Including these processes would drastically alter the physics of the model- 
essentially, from the quantum mechanics of one variable to quantum field theory— 
and would require a formalism for describing the topology changes. It should not, 
however, alter the instability suggested by our results. At issue is the stability 
of the spacetime foam, of which an individual wormhole is just one fluctuation. 
Certainly for any particular wormhole, the probability of growth to large size is 
affected strongly by the inclusion or exclusion of topology change. Indeed, since the 
time scale implied by our results for the wave function to leak to large radii is much 
longer that that for a classical bounce, the chance that a specific wormhole grows 
large should be much smaller than that it pinches off and disappears at zero radius, 
if that is allowed with more than an extremely small probability. But given the 
possibility of topology change, the foam should contain an equilibrium population 
of holes fluctuating into and out of existence. If it is possible for a hole to grow 
large, this population will eventually give rise to some large holes. Again the analogy 
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may be drawn to the a decay of a heavy nucleus: a particles continually form and 
disperse within the nucleus, on a time scale typically much shorter that that of 
the decay; the instability represented by the tunnelling of an a particle out of the 
nucleus remains. 

The more sophisticated analyses needed to probe the quantum dynamics of 
spacetime beyond the restrictions and limitations of these calculations present a 
considerable challenge. Our results suggest, however, that this is an aspect of 
quantum gravity theory well worth such consideration. 
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Figure Captions 

FIG. 1. Amplitude of the wormhole propagator G?( WKB ) f, R , 0], for #0 = 
10 and t = 20, all quantities in Planck units. 

FIG. 2. Squared magnitude of a wormhole wave function ij)(R,t) evolving via 
the propagator 

G (WKB)_ The initial 

wave function is a real Gaussian, as described 
in the text. All quantities are in Planck units. 

FIG. 3. Magnitude of the dominant contribution ip to a wormhole wave func- 
tion taken to be initially unity near the radius Rq = 10 Planck lengths, evaluated 
at the same radius after 10 4 -10 6 Planck times. 

FIG. 4. Behavior of the WKB propagator at large imaginary time r, here in 
units of the Planck time. The initial and final radii are fixed at 10 Planck lengths. 
Shown are the real and imaginary parts of the r derivative of In G. If the Feynman- 
Kac ground-state energy [Eq. (4.4) of the text] existed, the real part would approach 
that value and the imaginary part would tend to zero at large r. 
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